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Abstract 

We consider the elliptic system of linear elasticity with bounded mea- 
surable coefficients in a domain where the second Korn inequality holds. 
We construct heat kernel of the system subject to Dirichlet, Neumann, or 
mixed boundary condition under the assumption that weak solutions of the 
elliptic system are Holder continuous in the interior. Moreover, we show 
that if weak solutions of the mixed problem are Holder continuous up to 
the boundary, then the corresponding heat kernel has a Gaussian bound. In 
particular, if the domain is a two dimensional Lipschitz domain satisfying 
a corkscrew or non-tangential accessibility condition on the set where we 
specify Dirichlet boundary condition, then we show that the heat kernel 
has a Gaussian bound. As an application, we construct Green's function 
for elliptic mixed problem in such a domain. 
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1 Introduction 



In a domain (i.e. a connected open set) Q c R" (n > 2), we consider the 
differential operator 

where h is a column vector with components u 1 , . . . , u n , A a P(x) are n x n 
matrices whose elements cl^(x) are bounded measurable functions which 
satisfy the symmetry condition 



«#(x) = i*{x) = «J(x), (1.2) 
S T f <a%(x)de a <\K 2 \s + ?\ 2 , (1.3) 



where (5 is an arbitrary n x « matrix with real entries £,", Ki,k 2 > 0, and 
141 = {£,, £) 1/2 . Here and below, we will follow the convention that we 
sum over repeated indices and notation (£, rf) = £,' a rf a for n x n matrices 
E, = (E,' a ) and r\ = (rf a ). The operator L defined by <l.lt can also be written in 
coordinate form as follows. 

M-±l£;(4U%). <- 

a,p=l ;=1 " \ P ' 

Our assumptions on the coefficients a"f will include the equations of linear 
elasticity and in this case, the coefficients a°? (x) are usually referred to as the 
elasticity tensor; see e.g. 1 2 . 18 1 . In the classical theory of linear elasticity, the 
elasticity tensor for a homogeneous isotropic body is given by the formula 

a". = ASiadjp + li(5ij6 a p + bipbja), 

where A > 0, j.i > are the Lame constants, 5ij is the Kronecker symbol. 
In this case, the conditions l ll.2t and <1.3t are satisfied with K\ = 2j.i and 
k 2 = 2[i + n/\. We define the traction t = t(m) by the formula 

r(u) = v a A^(x)p-, (1.4) 

OXp 

where v = (vj, . . . , v n ) is the unit outward normal to <9Q. 

We consider following boundary value problems, which are the ones 
most frequently considered in the theory of linear elasticity. 

1. Dirichlet (displacement) problem 

Lu = / in Q, u = O on <?Q. (DP) 

2. Neumann (traction) problem 

Lu = / in Q, t(m) = cp on <?Q. (NP) 



2 



3. Mixed problem 



In = / in Q, u = O on D, t(u) = q> on N, 



(MP) 



where D and N are subsets of <9Q such that D U N = <9Q and D n N = 0. 

In the above, the equation as well as the boundary condition should be 
interpreted in a weak sense; see Section[2]for a precise formulation. 

In this article, we are concerned with the heat kernel associated with 
the mixed problem (MP). By allowing D = <9Q,N = and D = 0,N = dO. 
in JMP1 , we may regard dDFl and JNPt as extreme cases of (MP) and this 
is why we focus on (MP). By the heat kernel for (MP), we mean an n x « 
matrix valued function K(x, y, t) satisfying 



where 6 y (-) is Dirac delta function concentrated at y, I is the nxn identity 
matrix, and the equation as well as the boundary condition should be inter- 
preted in some weak sense. The heat kernels for (DPI and (NP) are similarly 
defined and they are frequently referred to as Dirichlet and Neumann heat 
kernels. 

We assume that Q is a bounded (e, <5)-domain of Jones fl5| and if D # 
and D <9Q, we assume further that it has a Lipschitz portion; if D = <9Q, 
then we require none of these conditions (see HI in Section |2~5) . We prove 
that if weak solutions of the system Lit = are locally Holder continuous 
(see H2 in Section |2~51 , then the heat kernel for (MP) exists and satisfies a 
natural growth estimate near the pole; see Theorem 13. II It is known that 
weak solutions of L are Holder continuous if n = 2 or if the coefficients are 
uniformly continuous. We also prove that if the gradient of weak solutions 
of the system Lu = satisfy the growth condition called Dirichlet property 
(see H3 in Section I2IB1 , then the heat kernel has a Gaussian upper bound; 
see Theorem 13.101 The Dirichlet property is known to hold in the case 
when Q is a Lipschitz domain in R 2 and D satisfies the corkscrew condition 
(see 1 25 1) or when the coefficients and domains are sufficiently smooth and 
D n N = (see [9]). As an application, we construct Green's function for the 
elliptic system from the heat kernel and in the presence of Dirichlet property, 
we show that the Green's function has the usual bound of C\x - y\ 2 ~" (or 
logarithmic bound if n = 2); see Theorem l4.3l 

A few remarks are in order. The Dirichlet or Neumann heat kernels for 
elliptic equations are studied by many authors; see Davies [6], Robinson 
(23), Varopoulos et al. l27l , and references therein. For the heat kernel 
for second-order elliptic operators in divergence form satisfying Robin- 
type boundary conditions, we mention Gesztesy et al. flOl . Dirichlet and 
Neumann Green's functions for strongly parabolic systems are studied in 
Cho et al. |3 4j and Choi and Kim \5\. We should mention that our paper, 
though technically more involved, is an extension of their method. The 
elliptic Green function for (MP) in two dimensional domains is constructed 
in Taylor et al. |25| by a different method not involving the heat kernel. 



§- t K(x,y,t)-L x K(x,y,t) = 
K(x, y,t) = 
r x (K(x,y,t)) = 
K(x,y,Q) = b v I 



in Q x (0,oo), 
on D x (0, oo), 
on N x (0, oo), 
on Q, 



(1.5) 



3 



Previously, Taylor et al. |26) give a construction of the Green function 
for a class of mixed problems for the Laplacian in a Lipschitz domain in 
dimensions two and higher. 

In recent years, there has been increasing interest in the study of elliptic 
equations under mixed boundary conditions from a variety of viewpoints. 
Haller-Dintelmann et al. fl4l consider Holder continuity of solutions of 
a single equation with an interest in applications in control. Mazzucato 
and Nistor [20] study the mixed problem for the elliptic system of elas- 
ticity in polyhedral domains under mixed boundary conditions. Finally, 
a recent monograph of Maz'ya and Rossmann fl9l treats mixed problems 
for systems in polyhedral domains. We refer to the above works and their 
references for background on the study of mixed problems for systems. It 
would be interesting to see if the well-posedness results in these works can 
be used to obtain fundamental estimates H2 and H3 that we use to obtain 
further estimates for the heat kernel. 

The organization of the paper is as follows. In Section[2] we introduce 
some notation and definitions including the precise definition of the heat 
kernel for (MP). In Section [3] we state our main theorems (Theorems 13.11 
and !3.10V which we briefly described above. In Section|4j we construct, as 
an application, the Green's function for jMPl and obtain the usual bounds. 
We give the proofs for our main results in Section [5] and some technical 
lemmas are proved in Appendix. 

2 Preliminaries 

2.1 Notation and definition 

Throughout the article, we let Q denote a domain in R" and let D and N be 
fixed subsets of <?Q such that DUN = dQ and D n N = 0. We use X = (x, t) 
to denote a point in R" +1 ; x = {x\, . . . ,x n ) will always be a point in R". We 
also write Y = (y, s), X = (x , to), and reserve notation 

Y=(y,0) = (y 1/ ...,y„,0). 

We define the parabolic distance in R" +1 by 

|X- Y\g* = max(|x - y\, j\t - s\), 

where | • | denotes the usual Euclidean norm, and write |X| a? = |X - 0| . 
We use the following notation for basic cylinders in R" +1 . 

Q(X,r) = \Y e R" +1 : |Y-X|^ < r\, 
Q_(X,r) = {Y=(y,s)eR" +1 : \Y-X\g* <r,s< t), 
Q + (X,r) = {Y = (y,s)eR" +1 : \Y-X\& <r,s>t}. 

We also use B(x,r) = [y e R" : \y - x\ < r\ to denote a ball in R". For the 
vector valued function u = (m 1 , . . . , n") T , we denote by e(u) (called the strain 
tensor) the matrix whose elements are 

, , 1 1 die dui\ 

e <; (M)= 2fe + ^-J- 
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We denote by ft the linear space of rigid displacements of R"; i.e. 

ft = {Ax + b : A is real n x n skew-symmetric matrix and b e R"l (2.1) 
Note that ft is a real vector space of dimension N = n(n + l)/2 and that 

(it, v) = I u ■ v dx = I v T u dx (2.2) 
Jo Jq 

defines an inner product on ft. Let us fix an orthonormal basis {(Oi}f =1 in ft 
and define the projection operator tz« : L 2 (Q)" — > ft by 



(«) = £<*, <w,-><Wi. (2.3) 



7T« 



The above formula still makes sense for u = b y ek, where e/ c is the fcth unit 
(column) vector in R". For y e Q, we denote by T y = T y (x) an n x n matrix 
valued function such that 

T ti»M ifD = 0, 
y |0 if D * 0. 

Roughly speaking, T y is an orthogonal projection of 5 f/ I on ^ if D = and 
otherwise. We set 

a s 

B3(ii,v) = a..- — - — . 

') dxp dx a 

It follows from <1.2t the form is symmetric (i.e., SS(u, v) = SS(v, it)) and 
from O) that 

k?38{u,u) < \e{u)\ 2 < x^SS{u,u). (2.5) 
It is easy to verify that O, O} imply (see |22] Lemma 3.1, p. 30]) 

^^n' a <\^ + ^\\n + n T \<^H- 

Finally, we denote d x = d(x) = dist(x, Q c ) when Q is clear from the context 
and write a A b = min(fl, b) and aV b = max(«, b) for a, b e R. 

2.2 Function spaces 

We use the notation in Gilbarg & Trudinger 1 13 1 for the standard functions 
spaces defined onUc R" such as U(U), W k 'F(U), C k ' a {U), etc. For I c dU, 
let W 1 ' 2 (Lf;r) be the subspace obtained by taking the closure in W 1,2 ((I) of 
smooth functions in U which vanish in a neighborhood of Y. Note that we 
have \N l;l {U;dU) = Wl' z (U). We shall denote 



/1,2 



w 1 ' 2 (u ; r) ifr# 

[it e W 1 ' 2 (Li) : f u udx = 0) if T = 



For Q and D as above, we define 



fw^CW ifD#0 
'{it € W^P": <u,z>> = 0, Vveft] if D = 0, ( ' ' 
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where ft and (u, v) are as in (2j) and (Z2j. Notice that V c W 1 ' 2 (Q; D)" . For 
fi € (0,1], we denote 

\u(x) - u(y)\ 

Wj&L = [u\, tU + |M|o;U = SUp — — + SUp \U{X)\. 

x,yeU \ X 3/1 xeU 

For spaces of functions defined on Q c R" +1 , we borrow notation mainly 
from Ladyzhenskaya et al. 1171 . To avoid confusion, spaces of functions 
defined on Q c R' !+1 shall be always written in script letters. For cj>\, we 
let Jf q (Q) denote the Banach space consisting of measurable functions on Q 
that are ^-integrable. For Q = Q x (a, V), we denote by «Sf 9 ,.(Q) the Banach 
space consisting of all measurable functions on Q with a finite norm 



I'd 1 ""' 



>•/</ 

t)\idx\ dt 



where q > 1 and r > 1. Thus JSf M (Q) is the space J?,(Q). By V^ IZ {Q) 
we denote the set of all bounded measurable functions u on Q for which 
|M|(i,fi/2 ; Q is finite, where we define the parabolic Holder norm as follows: 

\u\fi /f i/2;Q = [M](i,fi/2;Q + Mo;Q 

= sup K f-"; y)l + sup|u(X)|, fi E(0,l]. 

X*Y 

We write m € %°°(Q) (resp. ^ C °°(Q)) if U is an infinitely differentiable function 
on R" +1 with compact support in Q (resp. Q). We write D,-m = du/dxi 
(i = 1, . . . , n) and u f = du/dt. We also write Dn = D x u = (Diu, . . . , D„w). We 
write Q(fo) for the set of all points (x, t ) in Q and 7(Q) for the set of all t such 
that Q(t) is nonempty. We denote 

|m|q = I |D Y w| 2 dxdf + esssup I |x((x,f)| 2 dx. 

The space (Q) denotes the Banach space consisting of functions u 6 
-^(Q) with weak derivatives D,w e J£ q (Q) (i = 1, . . . , n) with the norm 

ll«ll^i.o (Q) = ll"IU,(Q) + I|D. v m||^( Q) 
and by W q (Q) the Banach space with the norm 

ll"ll^-U (Q) = NI.S?,(Q) + l|D.v«ll^,(Q) + l|WtlLsf,((3)- 

In the case when Q has a finite height (i.e., Q c R" x (-T, T) for some T < oo), 
we define %{Q) as the Banach space consisting of all elements of ^ 1,0 (Q) 
having a finite norm ||m||-^ 2 (Q) := \u\q and the space ^ 1,0 (Q) is obtained by 
completing the set W~ (Q) in the norm of ^(Q). When Q does not have 
a finite height, we say that u e ^(Q) (resp. ^ 1,0 (Q)) if U e y 2 (Qr) (resp. 
^ 1,0 (Qt)) for all T > 0, where Q T = Qn{\t\< T\, and |u| Q < oo. Note that this 
definition allows that 1 e i^ fi (Q x (0,co)). Finally, we write u e £? q ,ioc(Q) if 
u e J?,(Q') for all Q' <= Q and similarly define ^'"(Q), etc. 
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2.3 Weak solutions 

For /, g a e L 2 (U) n , where a = l,...,n, we say that u is a weak solution of 
Lu = f + D a g a in U if u e W^fli)" and for any v e W^U)" satisfies the 
identity 

[ ®{u,v) = - f f-v+ f g a -D a t>. (2.7) 
Ju Ju Ju 

Let r D , r N be disjoint subsets of dU. Recall that the traction t(m) is defined 

by the formula <1.4t . We say that m is a weak solution of 

Lu = f + D a g a in U, u = on T D , r(u) = g a V a on r N 

if u € VV^CHTd)" and for any z; e W^fU^U \ r N )" satisfies the identity 
<2.7> . Let Q,D,N be as above. For / e L 2 (Q)", we say that u is a weak 
solution of the mixed problem 

Lu = f in Q, « = on D, t(u) = on N 

if u e W 1 - 2 (Q; D)" and satisfies for any v e W 1 - 2 (Q; D)" the identity 

I ^(h,z>) = - I f-v. 
Ja Ja 

Let Q be a cylinder U x (a, fo), where U c R" and — oo < a < b < oo. For 
/ € £?2,i (QT and g a e j2f 2 (Q)", we say that u is a weak solution of 

"I -Lu = f + D a g a 

if u e %(Q)" and satisfies for all e ^"(Q)" the identity 

- f r s${u,$)= r r ^-d^. (2.8) 

Jq Jq Jq Jq 

Let Yd, Tn be disjoint subsets of dU. We say that u is a weak solution of 

u, -Lu = f + D a g a in Lf x (a, b) =: Q 
u = on r D x (a, b) 

= -#„v a on r N x (a,b)=:S 

if 1/ 6 %{Q) n , u(-,t) e W 1 ' 2 (!i;r D )' ! for a.e. t e (a,&), and it satisfies the 
identity <2T8j for all <b e ^™(Q U S) n . Next, denote Q = Q x (a,b), and let 
/ e j£f 2 ,i (Q)", g a e JS? 2 (Q)", and i/> e L 2 (Q)". By a weak solution in ^ 2 (Q)" 
(resp. ^'"(Q)") of the problem 

u t - Lu = f + D„g a in Q x (a, b) 

u = on Dx(a,fo) 

T ( M ) = -£« v « on Nx (a,b) 

w(-,a) = t/' on Q, 

we mean u(x, t) in f z {Qf (resp. ^'"(Q)") such that u(; t) e W 1 ' 2 (Q; D)" for 
a.e. f € (rt, b) and satisfying the identity 

- I «■</>,+ I ®{u,<\>)- [/</> + f g a -D a (f)= I tp(x)-4>(x,a)dx 
Jq Jq Jq Jq Ja 
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for any (f)(x,f) e ^.°°(Q)' ! that vanishes on D x (a,b) and equals to zero for 
t = a. We may identify / e .S?2,i(Q)" as an element in t}{a, b; V) in the sense 



lf(;t)](v)= f f(x,t)-v(x)dx, 



veV, 



J a 



and consider the problem 



Ut — Lu = f in Q x (a, b) 

u = on D x (a, b) 

t(m) = on N x (a, b) 

u{-,a) = xp on Q 



(2.10) 



u(; t) e V for a.e. t e (a, b). 



In the above, we impose the compatibility condition for ip that (xp,v) = 
for all v e ft in the case when D = 0. We shall say that m is a weak solution 
of the problem j2.101 if u 6 ^ 1,0 (Q)", u(-, f) e V for a.e. f e (a, b), and satisfies 
the identity 



for any t) e 'tf ™(Q) H that vanishes on D x (a, b), satisfies ($(•, t), »> = 
for all t e [a, b] and any v e % and equals to zero for t = a. Note that if 
D 0, then a weak solution of the problem <2.10t is also a weak solution 
in ^ 1,0 (Q) B of the problem (Z9j with g a = for all a = 1, ... ,n, and vice 
versa. However, when D = 0, they are not the same in general; note that if 
/("/ _ /("/ e ^/ then they are the same as elements in L 1 ^, b; V). 

2.4 Heat kernel for the system of linear elasticity 

We say that an n x n matrix valued function K(x, y, f), with measurable 
entries Kij : Q x Q x [0, oo) — » R, is the heat kernel for l IMPt if it satisfies the 
following properties, where we denote Q = Q x [0, oo). 

a) For all y e Q, elements of K(; y, ■) belong to n ^(Q \ Q+(f , r)) 
for any r > 0. 

b) For all y e Q, K(-, y, •) is a generalized solution of the problem <1.51 
in the sense that K(-,y,t) e W 1 ' 2 (Q;D)" for a.e. t > and for any 
4> = (cp\ ...,(p") T e %°°(Q)" that vanishes on N x [0, oo), we have the 
identity 



Jq Jq Jq Jq 



i/<(x) ■ <p{x, a) 



dx 




+ a"j s —K jk (x,y,t)—<p i (x,t)dxdt = (j> k (y,0). (2.11) 




c) For any / 



(/\ . . . ,f"f e tf™(Q)", the function u given by 



w 




Jo Jo 
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is, for any T > 0, a unique weak solution in "/^ ' (Q x (0, T))" of the 
problem 

w f -Lu = f in Q x (0, T), 
u = on D x (0, T), 

t(m) = on N x (0, T), 

m(-,0) = on Q. 



(2.12) 



We note that part c) of the above definition gives the uniqueness of the heat 
kernel for <MP|. 



2.5 Basic Assumptions and their consequences 

HI. We assume Q c R", n > 2 is a bounded domain. If D = <?Q, we do 
not make any further assumption. Otherwise, we assume the (e, 5)- 
condition of Jones fl5l for some e, 5 > 0: For any x, y e Q such that 
|x - y\ < b, there is a rectifiable arc y joining x to y and satisfying 

1 e\x — z\\y — x\ 
l(y) < -\x — y\, d(z) > j j for all z on y, 

^ \% ~ y\ 

where l(y) denotes the arc length of y and d(z) is the distance from z 
to the complement of Q. If D # <?Q and D + 0, we assume further 
that D has a portion of Lipschitz boundary; i.e. there exist x e D 
and a neighborhood V of x in R" and new orthogonal coordinates 
{yi, . .. ,y„] such that V is a hypercube in the new coordinates: 

V = 1(1/!, . . . , y„) : -a, < xjj < a,, l<j< n); 

there exists a Lipschitz continuous function <p defined in 

V = {(yi, . . . , y n -i) : -Qj < xjj < ay, 1 < ;' < n - 1}; 

and such that 

|<p(y')l < «n/2 for every y' = (y u y„-i) e V , 
O n V = \y = (y', y») eV:y„< <p(y')l, 
Dny=jy = (y',y„)ey:y„=<p(y')). 

In other words, in a neighborhood V of x , Q is below the graph of (p 
and D is the graph of cp. 

Basically, we introduce the assumption HI is to guarantee the multiplicative 
inequality J2.15t and the second Korn inequality J2.18t are available to us. 
We recall that the following multiplicative inequality holds for any u in 
W!' 2 (R") with n > 1; see El Theorem 2.2, p. 62]. 

INb^p^ < c(n)||D M [|«|| M [|W. (2.13) 

If we assume HI, then there is an extension operator £ : W 1,2 (Q) — > W 1/2 (R") 
such that the following holds; see l24l Theorem 8]. 

||£m|| L 2 (r „) < C||M|| L 2 (n) , HEullwUdj,,) < C||w|| (2.14) 
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Then by combining ( 12.141 and J2..13I , for any u e W 1,2 (Q; D), we obtain 

\\u\\„ {D) < C||D(E M )||^^||E«|g;^ 

where in the last step we used HI to apply the Friedrichs inequality (or 
Poincare's inequality if D = 0): 

Nbto, < C||DM|b (D) , V«eW u (Q;D). 

We have proved that HI implies that there is y = y(n, Q, D) such that for 
any u e W 1 ' 2 (Q; D), we have 

IMI^CT^yllDull^lN^. (2.15) 

If u e r 2 (0 x is such that u(-,f) e W 1,2 (Q;D) for a.e. t e («,&), where 
-co < a < & < oo, then by 1 12.151 we have (see [jZl pp. 74-75]) 

l|w||^P2<„ + 2)/»<« x ("») - yNoxW)- (2.16) 

Another important consequence of the inequality (2.151 is the following: 
For any u e ^z(Q x (a, b)) with b - a < oo, we have 

IMI^ 2(n+2)/ „(ax W) ) < {2y + Q>- fl)"' 2 ^!- 1 )^ | m |qx W) . (2.17) 

We refer to [17] Eq. (3.8), p. 77] for the proof of 12.171 . Moreover, HI implies 
the following second Korn inequality: (see (8) for the proof) 

||«|| W U( Q) < C {\\u\\ LHa) + IkWIbjq)} . (2.18) 

In fact, if u e W 1 - 2 (Q;<9Q)" = W,5' 2 (Q)", we have the first Korn inequality 

l|DM|| 2 2(n) <2|| £ („)|| 2 2(Q) . 
Also, we have the following inequalities for any u e V: 

I|h|| wU{q) < C||e(u)|| L2(Q) . (2.19) 

The inequality (2.191 is obtained by utilizing {2.181 in the proof of l22l 
Theorem 2.7, p. 21]. By (2.191 and (23), for any w 6 V, we have 

I m{u,u)dx>c I \Du\ 2 dx. (2.20) 

Lemma 2.21. Assume HI and let ip 6 L 2 (Q)' ! and f e «Jz? 2 ,i (Q)", Q = 

Q x (a, b) and -oo < a < b < oo. Twen, i/zere exists a unique weak solution u 
in ^ 1,C (Q)" of the problem (2.91 . Moreover, if we assume that (\p,v) = for all 
v eft in the case when D = 0, then there also exists a unique weak solution of the 
problem (2.101 . l/ll/ll.sf 2( „ +2)/( „ +4) (Q) < °°, then the weak solution u of the problem 
(2.101 satisfies an energy inequality 

+ \\rp\\ LHa ^, (2.22) 
where C depends only on n, Kj, jc 2 and the constants appearing in (2.161 and (2.191 . 
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Proof. With the aid of the second Korn inequalities 42.18t or J2.191 , it follows 
from the standard Galerkin's method and the energy inequality. ■ 

H2. There exist p e (0, 1] and A > such that if u is a weak solution of 
Lu = in B = B(x , r), where x e D. and < r < d{x ), then u is Holder 
continuous in jB = B(x ,r/2) with an estimate 

[u] m: t B <A r-^U\u(y)\ 2 dy\ . (2.23) 

Here, we use the notation ■£ u dx = ^ J g udx. 

It follows from H2 that a weak solution of u t - Lu is also locally Holder 
continuous. 

Lemma 2.24. H2 implies that there exist e (0, fz ) an d Ai > such that 
whenever u is a weak solution in %(Q)" ofu t - Lu = in Q = Q-(X ,r), where 
Xo = (x , t ) e Q and < r < d(x ), u is Holder continuous in \Q = Q_(X , r/2) 
and we have an estimate 

Proof. With the second Korn inequality available to us, the proof is essen- 
tially the same as that of 1 16 Theorem 3.3]. ■ 

Finally, we introduce a condition that was originally considered by 
Auscher and Tchamitchian [l] and is referred to as the Dirichlet property. 

H3. There exist p e (0, 1] and A > such that if u is a weak solution of 

( Lu = in BnQ, 
u = on BnD, 
[ t(m) = on BnN, 

where B = B(x , r) with x e Q and < r < diam Q, then for any 
< p < r, we have 

( \Du\ 2 dx< A (-) '° f \Du\ 2 dx. (2.25) 

The following lemma says that H3 implies H2. Moreover, it shows that if 
there is p > such that for all x e Q and < r < diam Q, we have 

|Q n B(x , r)\ > fir' 1 , (2.26) 

then, weak solutions of w, - Lu = with homogeneous boundary data are 
Holder continuous up to the boundary. It is not hard to check that (e, 6)- 
domains satisfy condition 12.26t , so that domains that satisfy HI also satisfy 

em. 

Lemma 2.27. Let Q = Q x [0,oo), D = Dx [0,oo), and N = N x [0,oo). If 
Q satisfies the condition H.261 , then H3 implies that there exist e (0, p. ) and 
A\ > such that ifu is a weak solution in f^iQ n Q)" of 

{ ti t -Lu = in QnQ, 

J m = on Q n £), (2.28) 

I t(h) = on QnM, 
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where Q = Q-(X , r) with X = (x , t ) eQ and < r < yftg A diamQ, then u is 
Holder continuous in |Q n Q, where \Q = Q-{X , r/2), and we have an estimate 

r^M^iQna + l"lo ;iQ n^ * A 1 r-^\\u\\x 2(QnQ) . (2.29) 
Proo/. See Appendix l6.ll ■ 

3 Main theorems 

Theorem 3.1. Assume the conditions HI and H2. Then there exists a unique 
heat kernel K(x, y, t)for <MPt . It satisfies the symmetry relation 

K(x,y,t)=K(y,x,t) T (3.2) 

and thus by J2.12> , for any f e ^"(Q x [0, £»))", we have 

u(x,t)= I I K(x,y,t - s)f(y,s)dyds (3.3) 
Jo Jd 

is, for any T > 0,a unique weak solution in % ' (Q x (0, T))" of the problem 42.12t . 
Also, for ip € L 2 (Q)' ! , the function u given by 

u{x,t) = J K(x,y,t)il>(y)dy (3.4) 



is, /or flwy T > 0, a unique weak solution in ^ 1,0 (Q x (0, T))" o/f/ze problem 

(3.5) 



m, - Lm = in Q x (0, T), 

h = on D x (0, T), 

t(u) = on NX (0, T), 

w(-, 0) = i/> on Q 



And ;/i/> is continuous at x e Q in addition, then 

, hm m I K(x,y,t)ip(y)dy = tp(x ). (3.6) 



xeQ,(>0 



Moreover, the following estimates holds for all y e Q, zrfiere roe use notation 
Q = Q x [0, oo ), i y = %), and Y = (y, 0). 

J) ll^y,-)ll^ (Q+( f,r)) ^ C p r-'^'P, Vr e (0,d y ] ; Vp e [l, *f). 

2) \\(x,t) eQ: \K(x,y,i)\ > A}\ < CA- {n+1)!n , VA > d~ n . 

3) l|DxK(-,y,OII^(Q + (f / ))^C p r— 1+ C +2 )/'' / Vre(0,d y ] ; Vpe[l,{$). 

4) |{(x,f) € (3: IDxKfoy,^ > A}| < CAT**, VA > d""" 1 . 

5) For X = (x, t) e Q satisfying \X-Y\a? < d v /2, we have 

\K(x,y,t)\<C\X-Y\-£. (3.7) 
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6) For X = (x, t) and X' = (x', f) in Q satisfying 

2\X'-X\&<\X-f\&<d y l2, 

we have 

\K(x', y,t')- K(x, y,t)\< C\X' - X\% |X - Yf^' 1 , (3.8) 

In the above, C are constants depending only on n,Ki,K 2 ,Q,D, j.i ,A and C p 
depend on p in addition. 

Remark 3.9. It will be clear from the proof that besides the estimates 1) - 6) 
in Theorem l3.ll we also have 

7) WK{; V, OlUwoufcto) * Cr-n Vr e (0, d y \. 

8) \K(, y, -)\ QXQ4tr) < Cr-"' 2 , Vr e (0, d v ]. 
Here, we use the notation 



K(x, y, t) 



K(x, y,t) if D 1 i 

K(x,y,t)-n K {b y I){x) if D = i 



where Tin is as defined in J2.3I ; see also J5.28L Moreover, if Q is such that it 
admits a bounded linear trace operator from W 1,2 (Q) to L 2 (<9Q), then it can 
be shown that for / e (Q x (0, T))" and g e jS* (N x (0, T)f, where <?* 
and rk (k = 1,2) are subject to the conditions of [17| Theorem 5.1, p. 170], 
the function u defined by 

u(x, t) = I I K(x, y, f - s)/(y, s)dyds+ I I K(x, y, t - s)g(y, s) dS y ds 
JoJa JoJa 

is a unique weak solution in f^' (O x (0, T))" of the problem 

' J*ftt = / in Q x (0, T) 

w = on D x (0, T) 

t(m) = g on N x (0, T) 

u(-,0) = on Q. 

The proof is similar to that of the representation formula {33} and is omitted. 
Theorem 3.10. Assume HI and H3. There exists the heat kernel K(x, y, t) for 
{MP} and it satisfies all the properties stated in Theorem \3J] Moreover, for x, y e Q 
and t > 0, we have the Gaussian hound 

C (-9\x-y\ 2 ) 

\K{x,y,t)\<— -exp ' Jl , (3.11) 

(VfAdiamQ) I 1 I 

where C = C(n,Ki,K 2 ,Q,D, j.i ,A ) and 9 = &(ki,k 2 ) > 0. Furthermore, for 
X = (x, t) and X' = (x' , f) in Q satisfying |X' - X| & < \ (|X - Y| & A diam o), 
we have 



\K(x',y,t')-K(x,y,t)\ < C 



where is as in Lemma \T27\ 



\X' —X\g? 



\X - Y| & A diam Q 
1 

( ^t A diam q) 



1 ( $\x - y\ 2 , 

x 7^- — ex P ( 3 - 12 ) 
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4 Applications 

4.1 Some examples 

1. If the coefficients are constant, then it is well known that H2 holds 
with j.i = 1 and A = A (n,K\,K 2 ). In fact, it is also known that if 
the coefficients belong to the VMO class, then H2 holds with jUo and 
A depending on the BMO modulus of the coefficients as well as on 
n, Ki, k 2 . Therefore, the conclusions of Theorem l3.1l are valid in these 
cases. 

2. If the coefficients belong to the VMO class, the domain Q is of class 
C 1 , and D n N = 0, then it is known that H3 holds. Therefore, the 
conclusions of Theorem |3.10| are valid in this case. 

3. If n = 2, then it is well known that H2 holds with j.i = /.<o(ki,K2) 
and A = A (k 1 ,k 2 ). Therefore, the conclusions of Theorem 13.11 are 
valid. In fact, if Q is a Lipschitz domain and D is a (possibly empty) 
set satisfying the corkscrew condition, i.e., for each x e dD (where the 
boundary is taken with respect to <?Q) and r e (0, r ), we may find 
x r € D so that \x - x r \ < r and dist(x r , <?Q \ D) > M"V, where r > 
and M > are constants, then it is known that H3 holds; see |25|. 
Therefore, the conclusions of Theorem |3.10| are valid in this case. 

4.2 Green's function for the elliptic system 

We say that an n x n matrix valued function G(x, y) is the Green's function 
of L for IMP) if it satisfies the following properties: 

i) G{; y) e Wj£(Q) and G(-, y) e \ B(y,r}) for all y e Q and r > 0. 
In the case when D = 0, we require J n v(x) T G(x, y)dx = for any 
veR. 

ii) G(-, y) is a weak solution of 

-LG(-,y) = 6yl - T y in Q, G(-,y) = on D, r(G(; y)) = on N 
in the sense that we have the identity 



for any (p = (cp 1 , . . .,</>") r e C~(Q)" n V; see (H3) and for the 
definition of T y and V. 

iii) For any / e C™(Q)" n V, the function u defined by 





(4.1) 



is the weak solution in V of the problem 




on 



m 



on 



Q, 
D, 
N. 



(4.2) 



14 



We note that with aid of the second Korn inequality d2.19t , which is valid 
for any u e V, the unique solvability of the problem <4.2t in the space V is 
an immediate consequence of Lax-Milgram lemma. We also note that the 
property iii) of the above definition together with the requirement 

J v(x) T G(x, y) dx = 0, VveK 

gives the uniqueness of the Green's function. 

Theorem 4.3. Assume the conditions HI and H2. Then, there exists a unique 
Green's function G(x, y) for tjMV\ . We have 

G(y,x) = G(x,y) T (4.4) 

and thus by ( 14.11 , for any f e C°°(Q)" n V, we find 




G(x,y)f(y)dy 



is a unique weak solution in V of the problem J4.2t . If we assume H3 instead of 
H2, then for any x, y e Q, we have 

K*rt 1C {l + h (f=£)}. (4.5, 

ii) n>3 

\G(x,y)\<C\x-yt'\ (4.6) 

and moreover, for any x, y e Q satisfying \x- x'\ < \\x- y\, we have 

\G(x', y) - G(x, y)| < C\x' - x\^ \x - yt"'^ . (4.7) 

In the above, C is a constant depending on the prescribed parameters and diam Q 
and \i\ e (0, 1) is as in Lemma \2~27\ 

Corollary 4.8. Let n = 2 and assume HI. Then, there exists the Green's function 
for iMF\ that satisfies J4.4t . Moreover, if CI is a Lipschitz domain and D satisfies 
the corkscrew condition described in Section 14.11 then the estimates {45}, J4.7I 
hold, and pi and C are constants determined by X\, k 2 , Q, and D . 

Proof. Follows from Example 3 in Section l4~T1 and Theorem l4.3l ■ 

5 Proofs of main theorems 
5.1 Proof of Theorem [37T1 

In the proof, we denote by C a constant depending on the prescribed pa- 
rameters n, Ki, k 2 , jJ-o,Ao as well as on Q and D; if it depends also on some 
other parameters such as p, it will be written as C p , etc. 

We fix a O e C™(R") such that O is supported in B(0, 1), < <J> < 2, and 
L„ <I> = 1. Let y e Q be fixed but arbitrary. For e > 0, we define 

OV(x) = e-"0((x-y)/e) 
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and let v e = v c ,y,k be a unique weak solution in f^°{Q. x (0,T))" of the 
problem 

U t - Lu = in Q x (0, T), 
u = on Dx(0,T), 

t(m) = onNx(OJ), ( ' 

u{-, 0) = Oy^e* on Q, 

where et is the k-th unit column vector in R"; see Lemma T2.21I By the 
uniqueness, we find that v e does not depend on a particular choice of T and 
thus by setting v e (x, t) = for f < and letting T — > oo, we may assume 
the c f is defined on the entire Q x (-oo, oo). We define the mollified heat 
kernel K e (x, y, t) to be an n x n matrix valued function whose fc-th column is 
v e , y ,k(x, t); i.e., 

K e jk {x,y,t) = v> e (x,t) = v[ Ajk (x,t). 

For / e ^f°(Q x (-oo, oo))", fix a, b so that a < < b and supp /cQx(u, b). 
Let u be a weak solution in ~f}' {Q, x (a, b))" of the backward problem 



-m ( - Lu = 
u = 
t(w) = 
u(-,b) = 



/ 



in Q X (a, b) 
on D x (a, b) 
on N x (a, b) 
on Q. 



Then, it is easy to see that we have 

f ® 1Jie (x)u k (x,0)dx= f f K e ik (x,y,t)f i (x,t)dxdt. 
Jn JoJn 

Next, we define K f (x, y, f) by 



K e (x,y,t) 



K e (x,y,t) ifD*0, 
K e (x, y, f) - l [0/0o) (f) n K (0 ;/ , e I)(x) if D = 0, 



(5.2) 



(5.3) 



where nR(0 1/e I)(x) is an n x n matrix whose fc-th column is 7T K (®, //e e<;)(x). 
We set v e = to be the fc-th column of K e (-, y, •). It is easy to verify that 
v e {-, t) e V for a.e. t > 0. Therefore, for any T > 0, it is the weak solution of 
the problem (see Section |23l and Lemma l2.21t 



Mf - Lu = 
u = 
t(m) = 
"(v0) = ^ yjJc 
M (-,t) 6 y 



where we denote 



e,y,k 



in Q x (0, T) 
on D x (0, T) 
on N x (0, T) 
on Q 

for a.e. t e (0, T), 



if D # 



(5.4) 



^Oy /e e, c -7T K (Oy /e e /c ) ifD = C 
By the energy inequality, we get (see Lemma [2.21> 



(5.5) 
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Let / be a smooth function supported in Q+(X ,R) c Q x (-00,00). Fix 
b > t + R 2 and let ft be the weak solution of the backward problem 



-u f - Lu = f in Q x (a, b) 

u = on D x (a, b) 

t(k) = on N x (a, b) 

u(-,b) =0 on Q 

u(-, t) e V for a.e. t e (a, b). 



(5.6) 



The unique solvability of the above problem is similar to Lemma l2.21l and 
by setting u{x, t) = for t > b and letting a — > -00, we may again assume 
that ft is defined on Q x (-00, 00). Then, similar to 1 12.221 , we have 



l»| 



Cll/ll. 



2 , 2(„+2)/{„+4)(Q+(X ,R))- 



nx(-oo,») 

By using Holder's inequality, we derive from | |5.7| that 

I|m||^ 2 (q + (Xo,r)) ^ CR 3+ " /2 ||/||j 5f „(Q +( x ,R)). 
Then by Lemma |5T9l below and the above estimate, we obtain 

l«lo;Q + (X ,R/2) < CR 2 ||/|L Sf<xl( Q + (x 0/ R))- 



(5.7) 



(5.: 



Lemma 5.9. H2 implies that u is continuous in Q + (X ,R/2) and satisfies the 
estimate 

|«| .i Q < c(R-<" +2 » /2 ||u||^ 2(Q) +R 2 ||/||^(o)), (5.10) 

where we denote aQ = Q+(X , <xR). In fact, the same conclusion is true if 'ft is a 
weak solution in %fQ) n of -u t - Lu = f (or u t - Lu = f) with f e Jf^Q)". 



Proof. See Appendixy 

Note that similar to (15-21 , we have the identity 

f O y , f (x)w ,c (x,0)dx= f f K^{;y,-)fdX. 
Ja J-00 Ja 

If B(y, e) x {0| c Q + (X , R/2), then fxFi) together with f5M yields 



(5.11) 



JJo + 



Kl(;y,-)fdX 



^ \u\o;Q + 



(X a ,W) 



< CR 2 



>Q+(X ,R) 

Therefore, by duality, it follows that we have 

\\k*(;y,-)\\.^(Q + (X , R)) <CR 2 



II.S?«,(Q + (X ,R))- 



(5.12) 



provided < R < d y and B(y,e) x {0} c Q + (X 0r R/2). For X e Q such that 
< d := |X - Y\ & < d y /6, if we set r = d/3, X = (y, -2d 2 ), and R = 6d, then 
it is easy to see that for e < d/3, we have 

B(y,e) x {0) c Q + (X ,R/2), Q.(X,r) c Q + (X ,R), 

and also that i> e = f> e ,yM is a weak solution in ^((^-(X, r)) of u t - Lu = 0. 
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Lemma 5.13. H2 implies that for any p > 0,we have 

Mtto-pwa) < C p r-^ 2) >"\\u\\^Q Am . 

Proof. It follows from the estimate J5.101 in Lemma 15.91 together with a 
standard argument described in fl2l pp. 80-82]. ■ 

Note that by Lemma l5.13l and <5.12t , we obtain 

\v e (X)\ < Cr-"- 2 ||z) e ||^ l(Q _ (X/)) < Cr-'- 2 \\v e ||^ i( q + (x ,r)> < Cd~", 

That is, for X = (x, t) e Q satisfying < |X - Y\ g? < d ;/ /6, we have 

\k £ {x,y,t)\<C\X-Y\-^, Ve<±|X-Y|^. (5.14) 

Next, we claim that for < R < d y , we have 

l^(v^-)b\ Q+( y,R,<CR-" /2 , Ve>0. (5.15) 

To prove 1 15.151 , we only need to consider the case when R > 6e. Indeed, if 
R < 6e, then (|53) yields 

\K e i;yr)\ Q \ Q+{ U) * \^i;yr)W-^) < Ce-' 2 < CR-' 2 . 

Fix a cut-off function C e ^f°(Q(f , R)) such that 

C = l on Q + (Y,R/2), < C < 1, \D X Q < 4R" 1 , IC/I < 16R' 2 . (5.16) 



By using the second Korn inequality j2.20| and <5.141 , we derive from j5.4| 
that 

sup \ \{l-Qv £ {x,t)\ 2 dx + ff \D x ((l - Qv e )\ 2 dxdt 

< C J£ (\D X Q 2 + |(1 - QCfl) \v e \ 2 dxdt 

< CR' 2 \\ |X - Y| -g dX < CR~", (5.17) 

which implies the desired estimate ( 15.151 . In fact, we obtain from \5.Y7\ that 
|(1 - QK e (-, y, Ob < CR-" 12 , Ve > 0. (5.18) 
We claim that for < R < d y , we have 

!!*•(•' 3/' •Jllawcoufc*")) " CR "" /2 - (5 ' 19) 
Ind eed, s et Q m := Q x (R 2 , oo) and Q p) := (Q \ B(y, R)) x (0, _R 2 ) and note that 
by 12.161 and 15.151 we have 

ll^(v^-)ll^ (n+2)/ „«3 (1) ) < CyW(;y,-)\ Qm < CyR-> 2 
and similarly, by 1 12.171 and J5.181 , we have 

ll^(v^-)ll^ (n+2)/ „(«3 (2) ) < C(2y + l)^yR-" 2 . 
By combining the above two inequalities, we get J5.191 . 



18 



Lemma 5.20. For any yeQ and e > 0, we have 

\\X=(x,t)eQ:\K e (x,y,t)\> A}\<CA-^, VA > d~ n , (5.21) 
|{X = (x,t) e Q: \D x K e (x,y,t)\ > A}| < CA"^, VA > (5.22) 
Also, for any y e Q, < R < d y , and e > 0, H?e tace 

\\K e (; y, -)\\^ QA U)) * C p R- + ™, Vp e [1, *±2), (5.23) 
||DJH, y, Oll^t*,, < C,,R-"- 1+ ™, Vp 6 [1, ). (5.24) 

Proof. We derive d5723t and J5.24t , respectively, from <5.21t and (5j22), which 
in turn follow from J5.19t and J5.15t , respectively; see (3) Lemmas 3.3 and 
3.4]. ■ 

Lemma 5.25. Suppose {ut}^ is a sequence in y 2 (Q) n such that supju^lg < 
A < oo, then there exists u e f^iQ)" satisfying \u\q < A and a subsequence 
that converges to u weakly in W~ (Q, x (0, T)) n for any T > 0. Moreover, if each 
uic(-, t) e V for a.e. t e (0, oo), then we also have u(-, t) e V for a.e. t e (0, oo). 

Proof. See (3] Lemma A.l]. ■ 

The above two lemmas contain all the ingredients for the construction 
of a function K(-, y, •) such that for a sequence e fl tending to zero, we have 

K^(-,y,-)^k(;y,-) weakly in W^(Q + (¥,dy)f , 

(1 - QK e » (•, y, •) — (1 - QK(; y, •) weakly in W 2 lfi (Cl x (0, T))"\ (5.26) 

where 1 < q < a *?, C is as in J5.161 with R = Jy/2, and T > is arbitrary. It is 
routine to check that K(-, y, •) satisfies the same estimates as in Lemma l5.20l 
as well as J5.15t and 45.191 ; see Section 4.2]. Note that by Lemma I5T251 
we have K(-, t,y) e V for a.e. t > 0. We define K(x, y, t) by 

[k(x, y, t) if D * 0, 

K(x, y,t):=i. * ' 5.27 
\K(x, y, t) + l [0/Oo) (0 n K {6 y I){x) if D = 0. 

where nn(b y T) is an n x n matrix valued function whose fc-th column is 

N 
i=l 

where cu, = (a>j, ai") T e K and I <u, • w, dx = 6 i; -. (5.28) 

Then, it is easy to see that K(x, y, t) satisfies the estimates 1) - 4) in Theorem 
13.11 because K(x, y, t) satisfies all of them as we noted above, and H is a 
finite dimensional vector space so that all norms over K are equivalent. For 
example, to see the estimate 3) in the theorem holds, observe that 

ll D<u ill^P(Q + (tr)) ^ ^illwW(B(y,r))IQ+(*>)l^ Z 

< cr' !/2+( " +2)/ P < C(diamQ)" /2+1 r-"- 1+(,,+2)/ ''. 
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Since nR(b y ek) € "R, it is a weak solution of u t — Lu = 0. Therefore, by 
repeating the proof for 15.14> , we find K(x, y, t) satisfies the estimate i3.7i , 
while the estimate i3.8i is obtained from {3.7\ and Lemma \l. 241 We have 
thus shown that K(x, y, t) satisfies all the estimates 1) - 6) in Theorem l3.ll 

We now prove that K(x, y, t) satisfies all the properties stated in Sec- 
tion |Z4] so that it is indeed the heat kernel for (MP). First, note that the 
property a) is clear from 1), 3) in the theorem and 8) in Remark [3T91 To 
verify the property b), first note that by \5.3l together with 1 15.261 and \5.27\ , 
we have 

K^(;y,-)^K(;y,-) weakly in W q lp (Q + (?,d v )f , 
(1 - QIC" (•, y, •) — (1 - QK(; y, •) weakly in W 2 lfl (Cl x (0, T))"\ (5.29) 

for any T > 0. Next, suppose <j) = (cp 1 , . .., cp n ) T is supported in Q x (0, T) 
and note that by J5.lt we have 



By writing (j) = rjij) + (1 - r[)<p, where rj e V^°(Q{Y, d y )) satisfying r\ = 1 
on Q(Y,d y /2), and using <5.29t , and taking j.i — > oo in the above, we get 
the identity <2.11t : see (5) p. 1662] for the details. To verify the property 
c), let us denote f(x, t) = f(x,-t) and let m be a unique weak solution in 
^ 1,0 (Q x (-T, 0))" of the backward problem 



By letting T — > oo, we may assume that u is defined on Q x (-oo, 0). Then, 
similar to d5.2L for t > 0, we have 



We note H2 implies, similar to Lemma |5T9l that it is continuous in Qx(-oo, 0). 
By writing / = £/ + (1 - Qf and use ( l5!29l to get 



If we set u(x, t) = u(x,—t), then it becomes be a weak solution in >2(Q x 
(0, T))" of the problem f2~12l , and thus by the uniqueness the property c) 
is confirmed. Therefore, we have shown that K(x, y, t) is indeed the heat 
kernel for ((MP)- 

Now, we prove the identity | |3.2| . Let 




-u t — Lu = f 
u = 
t(m) = 
k(-,0) = 



in Q x (-T,0) 
on Dx (-T,0) 
on Nx(-T,0) 
on Q. 





KiV, x, s) = & w (y, s) = v[ Xxl (y, t-s) = K b u {y, x,t-s) 



(5.30) 
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and 

^(y,s) = 3>^(y,f-s), 

where v biX j and <t> Xib are as above. Observe that v bjXi i(y,s) is, for any -T < t, 
a unique weak solution in % ' (£1 x (-T, t))" of the problem 



-Z7 s -Li7 = in Qx(-T,f) 

Z7 = on Dx(-T,t) 

t(v) = on Nx (-T,t) 

v(-, t) = i> Xib ei on Q. 



(5.31) 



Then, similar to J5.21 , we have 



f *«(.,x,0)<V = f ^,(-,y,0^. (5.32) 
Jo Jn 

By repeating the proof of [3 Lemma 3.5], we obtain J3.21 from J5.321 as well 
as the the following representation of the mollified heat kernel: 

K e {x,y,t) = J K{z,x,tf% h£ (z)dz. 



In particular, by continuity of K(-, x, t) and 1 13.21 , we have 

\m\K e (x, y, t) = K(x, y, t). (5.33) 



Now, we turn to the proof of the formula 1 13.41 . Let u be the weak 
solution in ^(Q x (0, J)f of the problem l3~51 Le t X = (x, t) e Q x (0, T) 
be fixed but arbitrary and let v b = v biX j be as in {530). Then, it follows from 
the equations 1 13.51 and J5.311 that for sufficiently small 5, we have 

f ^{yW b {y,0)dy= ( u l (y,t)O brX (y)dy. 
Jn Jn 

Therefore, by using J5.301 , we obtain 

f K e f(y,x,t)f(y)dy= f U '(y,t)% liX (y)dy. (5.34) 
Jn Jo 

By 1 15.151 and ( 15.331 with x in place of y, we find by the dominated conver- 
gence theorem that 



Um \ ^i{y,x,tW{y)dy= \ K u (y,x, 
Jn Jn 



tW(y)dy. 

By Lemma I5T91 we find that u is continuous at X = (x, t). Therefore, by 
taking the limit /.i — > oo in J5.341 and using ( 13.21 we obtain ( 13.41 . 

Finally, let « be a weak solution ^'°(Q x (0,T))" of the problem i3~5l 
and <p be a Lipschitz function on Q satisfying |V</)| < K a.e. for some K > 0. 
Denote 

7(0:= f e 2 ^'\u(x,t)\ 2 dx. 
Jn 
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Then I'(t) satisfies for a.e. t > the differential inequality 



Y 



(t) = -2 \^(u,u) + 2^ a f™rL u i\dx 




dx 



(5.35) 



Then, by repeating the argument in |3 Section 4.4], we obtain the formula 



5.2 Proof of Theorem [3TT01 

By Lemma [2.27l and a remark preceding it, we observe that the conditions 
HI and H3 imply the condition H2 and the condition (LB) in @). Then, 
by Theorem 13.11 the heat kernel K(x, y, t) exists and by using ((535) and 
repeating the proof of [4 Theorem 3.1] with R max = diamQ, we obtain the 
Gaussian bound 43TTV Also, by E9t , for X = (x, t) and X' = (x' , t') in Q 
satisfying 

2|X' -X\g> <r<r := |X-f|^ A diamQ, 



\K(x', y, t') - K(x, y, t)\ < C|X' - X\% r^ n+2 ^\\K{; y, -^q^q) (5.36) 



Note that the estimate 43.11| > implies that for < s < (diam Q) 2 , we have 



From the above the estimate, we obtain the estimate 43.12t by repeating the 
proof of (5] Theorem 3.7]. More precisely, we consider the following three 
possible cases. 



If |X' - X\g» < r /8, then we take r = r /4 in (536} and use (|537) to get 

\K(x', y,t')- K(x, y,t)\< C|X - X\% , 
which implies <3.12t . If r /8 < \X' - X\g? < r /2, then we have 

\x' -y\< 3r /2 and r /2 < <r < diamQ 
and thus, by 13.111 we get 

\K(x', y, f) - K(x, y, f)| < \K(x', y, t')\ + \K(x, y, t)\ < Cr'", 
which also implies 43.12t . 



43.61 . The theorem is proved. 



we have 




(5.37) 



i) Case \x- y\< Vi < diam Q: In this case, we have 

r = Vt = \X-Y\& and \x-y\ 2 /t<l. 
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ii) Case Vf < \x — y|: In this case, r = \x - y\ < diamQ. Similar to (U 
Eq. (5.22)], for all (z, s) e Q.(X, r /2) n Q, we have 

\K(z, y,s)\< Cr"> 2 exp [~&\x- y\ 2 /4t} . (5.38) 

If |X' - X| &> <r /4:, then we take r = r /2 in (536) and use d538l to get 

\K(x',y,t')-K(x,y,t)\ < C\X' -X\%r-^r n/2 exj?{-d\x-xj\ 2 /4t\. 

If r /4 < |X' - X| iS3a < r ll, then use {3~Tfl and {538) to get 

|K(x', y, f) - JC(x, y, f)| < Cf"" /2 exp {-5|x - y| 2 /4f} . 

Therefore, we also obtain J3.121 in this case. 

iii) Case diamQ < rf: In this case r = d := diamQ, and the desired 
estimate (13.12b becomes 

|JC(x', y, f) - K(x, y, t)\ < C\X' - X\% cT""'' 1 exp {-9\x - y\ 2 /4t} . (5.39) 

Since t > d 2 , for all (z, s) e Q- (X, r /2) n (3, we have 

exp {-9|z - y| 2 /s} < e s/4 exp {-% - y| 2 /2f} . (5.40) 

If |X' - X|^ < r /4, then we take r = r /2 in J536t and use (5M[ to 
obtain d539l . If r /4 < |X' - X| ^ < r /2, then by 01T1 and l |5l0l 

\K(x', y, f) - K(x, y, t)\ < Cd~ n exp {-9\x - y\ 2 /2t] . 

Therefore, we also obtain J3.121 in this case. 
The theorem is proved. ■ 

5.3 Proof of Theorem [431 

Assuming HI and H2, we construct the Green's function G(x, y) for {MP) 
as follows. Note that (2.191 implies that there is a constant g such that for 
any u £ V, we have 

IMb(a) < e>lk(")lb ( o)- (5.41) 

By utilizing (5.4H and following the proof of E Lemma 3.2], we get that 
for x, y e Q with x y, we have 



|K(x, y, f)l dt < oo, 

where K(-, y, •) is as in the proof of Theorem |3.1| We then define 

G(x,y):= f K(x,y,t)dt. (5.42) 
Jo 

Then the symmetry relation J4.4t is an immediate consequence of 1 13.21 once 
we show that G(x, y) is the Green's function. We shall prove below that 
G(x, y) indeed enjoys the properties stated in Section l4~2l Denote 

-t 



so that we have 



R(x,y,t) := I K(x,y,s)ds 
Jo 



G(x, y) = limK(x, y, t). 



23 



Lemma 5.43. The following holds uniformly for all t > and y e O. 

i) \\K(; y, OlboM,)) < C p (d) + ?)df~ n , Vp e [l, f). 

ii) ||£( v y, f)|| l2 ( n+ 2,/„ (Q \ BM) < C(r 2 + , Vr e (0, d y ]. 

m \\dk(; y, t)\y m4y)) < c p (d 2 y + f)d; J 1 -" + " lp , Vp e [l, *±f ). 

W ||DK(-, y, Olli 2 (Q\ B(y/ )) < C(r 2 + ^r" 1 ""/ 2 , Vr 6 (0, d y ]. 

Proof. See Lemma 3.23]. ■ 

By the above lemma, elements G !; (x, y) of G(x, y) satisfy 

Gi,{; y) e W U (Q) and G, 7 (-, y) 6 W 1 ' 2 (Q \ B(y, r)) for any r > 0. 

Recall that columns of K(-, y, t) are members of V; see Lemma 15.25 1 There- 
fore, in the case when D = 0, for any v e% we have 

I u(x) T K(x / y, t)dx = 

JQ 



and thus, we also have 



X 



v(x) T G(x,y)dx = 0. 



We have shown that G(x, y) satisfies the property i) in Section l4~2l For the 
proof of the property ii) in Section l4~2l we refer to Section 3.2]. Finally, we 
show that the property iii) in Section l42l also holds. Let / e C~(Q)" n V and 
u be defined by the formula J4.lt . The integral l l4.ll is absolutely convergent 
by the property i) of section l42l Similarly, Lemma |5.43| implies that 

v(x,t) := J k{x,y,t)f(y)dy 

is well defined. Observe that 

v(x,t)= f f K(x,y,s)f(y)dyds= f f K(x,y,t - s)f(y)dyds. (5.44) 
JoJn JoJn 

Therefore, we have 

]imv(x,t)= l G{x,y)f(y)dy = u(x) (5.45) 
f -»~ Jo 



V t (X, t) ■ 



f K(x,y,t)f(y)dy. (5.46) 
Jn 



By J5.271 , the assumption that / e V, and j3~3l , we find from 1 15.441 that v is, 

for any T > 0, the weak solution in ^ a '°(Q x (0, T))" of the problem 

v,-Lv = f in Q x (0, T) 

» = onDx(OJ) 

t(v) =0 on N x (0, T) 

v(-,0) = on Q. 
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By a similar reasoning, J5.461 and the representation formula d3.4t implies 
that v t is, for any T > 0, the weak solution in ^ 1,0 (Q x (0, T))" of the problem 
(35\ with ijj = f. Then, we have (see Eq. (3.42)]) 



\\v t (; f)[b (Q) < CeT*^ f ||/|b (Q)/ Vt > 0. (5.47) 
Observe that by <5.44> , <5.46> , and the assumption that / e V, we have 

v(-, f)eV and v t {-, t)eV for a.e. t > 0. 
Also, note that for any <p = {(f) 1 , ... ,(p") T e V and for a.e. t > 0, we have 

f a f^T-(-> t )J- dx = f f'ftdx- f v\(-t)(p'dx. (5.48) 
Jq 7 cfrjs c« [V J Q J q 

Then, by setting (p = v(-, t) in ( 15.481 and using CHI}, for a.e. f > 0, we have 

\\e(v{;t))\\l 1(a) < C(||/|| L2(Q) + Ib f (-,0lli2(c,)lb(v0lli 2( c, 

< C\\f\\ ma) \\e{v (;t))\\m a) , 
where we have used J5.471 . Therefore, by J2.191 . for a.e. t > 0, we have 

lb (vt)llwW(Q) < C|l/lb(D)- 

Then, by the weak compactness of the space W 1/2 (Q)' ! together with the fact 
that V is weakly closed in W 1/2 (Q)", we find that there is a sequence U»i)™ =1 
tending to infinity and u £ V such that 

v(; t m ) u weakly in W lz (Q) n . 

By J5.45b . we must have u = ueV and thus, for all (p e V, we get 

,. f apdvl , f afi dui 

lim fl.r- — (vfmH — ax = «.. - — - — dx. (5.49) 
Then, by <5T49> , <5T47t , and 0!48l , for any e V, we obtain 

r «°f dx = f f'cb' dx, 

J a '> dxgdx a J n J ' 



which shows u is a weak solution in V of the problem (42); see the remark 
that appears above Theorem |4.3l Therefore, we verified that G{x, y) defined 
by the formula J5.421 also satisfies the property iii) in Section l4~2l and thus 
it is indeed the Green's function for ((MP). 

Next, we assume H3 instead of H2 and proceed to prove the second 
part of the theorem. In the rest of the proof, we shall denote 

d := diamQ. 

By Theorem 13.101 we have the Gaussian bound 13.111 . In particular, for 
X = (x, f) e Q satisfying VI < diam Q, we have 

\K(x,y,t)\<C\X-Y\-£. (5.50) 
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Similar to (3) Eq. (6.17)], we have 

\k(x,y,t)\ < Cr-"e-^^ 2it - 2rl \ t>2r 2 , < r < d. (5.51) 
We set r := | mm(g, d). If < \x — y\ < r, then by ( |5.42t , we have 

/Hx-y| 2 p2r 2 poo 

\G{x,y)\< + + |£(x,y,i)|df =:Ii+I 2 +I 3 . (5.52) 

Jo J\x-y\ 2 J2r 2 

It then follows from j530t and <5.51t that 



7i < C f |x - y\~" dt < C\x - y\ 2 ~", 
Jo 

h<c f r" ,2 dt<\ 

Jl-v-J/l 2 { 



C + Cln(r/\x-y\) ifn = 2, 
C|x - y\ 2 -" if n > 3. 



h<C f r-"e- K ^ 2it - 2 ' 2) dt<Cg 2 r- n . 
Combining all together we get that if < \x — y\ < r, then 



|G(x,y)|<^ (5.53) 
1 C(l + (p/Jr) |x - y\ z " ifn>3. 



In the case when \x — y\> r, we estimate by <5.50> and <5.5H that 
|G(x,y)|< I \K(x,y,t)\dt + I |K(x,y,0ld£ 

Jo J 2 r 2 

2) 2 oo 

<C f r'"dt + C f r-"e- K ^ 2( '- 2rZ) dt<Cr 2 -" + Cg 2 r". (5.54) 

Jo J 2 r 2 



By J5.53b and <5.54t , we get l !4.5t and l !4.6t . Finally, we turn to the proof of 
the estimate ii.7) . Because we assume H3, the conclusions of Theorem l3.10l 
are valid. By <3l2l and the definition 0x27) , if |X - Y\ & < d, then we have 



\K{x', y, t) - K{x, y, t)\ < \K(x', y, t) - K(x, y, t)\ + C\x - x'\ 

< C\x' - x\ m \X - y|^" w whenever \x -x'\<\\x- y\. (5.55) 

We claim that for < r < d and t > 3r 2 , we have 

\K{x', y, t) - K(x, y, t)\ < C\x' - ^y-n-^^e-Ht-^) (5 56) 

whenever \x-x'\ < j\x-y\. Assume the claim J5.56t for the moment. Similar 
to J5.52t , in the case when < \x - y\ < r := \ min(p, d), we get 

\G(x',y)-G(x,y)\< 

Wx-yl 2 ^i 2 poo 

I + + \k(x',y,t)-k(x,y,t)\dt =:h+l 2 + h. 

Jo J|.v-)/l 2 J3r 2 
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It follows from <5.551 that 

Hx-yl 



rv-yr 

h<C\x'-x\^ \x - y\~"-^ dt < C\x' - x\^\x - y\ 

Jo 

h<C\x'-x\^ f r n/2 -^ l2 dt < C\x' - x\^\x - y\ 2 ~ n - 



fi 



Also, by J5.561 , we obtain 



Z 3 < Cr""-'' 1 |x' - x|* f e - K i ?" 2 <'- 2 ' 2 ) df 

< QV"-'' 1 |x' - x| w < C (| j 2 |x' - x|« |x - y| 2- "" 



Combining the above estimates together, we obtain 1471 when < \x-y\ < r. 
In the case when |x — y\ > r, by using ( 15.551 and J5.561 , we estimate 

\G{x',y)-G(x,y)\< f + f \K(x' ,y,t)-K(x,y,t)\dt 

JO Jir 2 

< C\x' - x^r 2 -"-^ + Cg 2 r"-^\x' - x| w 

< C(l + {Qlrf)\x' - x\^(rldf-"-^d 2 - n -^ 

< C(l + (g/r) 2 )(r/d) 2 -"-^\x' - x| w |x - j/! 2 ""-'' 1 . 



Therefore, we also obtain (47) when |x - y\ > r. It only remains for us to 
prove the claim ( 15.561 . The strategy is similar to the proof of ( 13.121 . Note 
that each column of K(-,y,t) is a weak solution in V 2 {Q) of u t - Lu = 
provided that Q <= Q \ |Y]. Therefore, similar to ( 15.361 , for < r < d and 
t > 3V 2 , we have 

\K(x', y, t) - K(x, y, f)| < C\x' - x|« r^-^m, V, OllaHQ-OWnfl) 
whenever |x - x'| < r/2. Then by ( 15.511 , we obtain j5.561 . ■ 

6 Appendix 

6.1 Proof of Lemma |2T27l 

We first show H3 implies H2. Suppose u is a weak solution of Lu = in 
B(z , r) c Q. By a well-known theorem of Morrey [21] Theorem 3.5.2], we 
have 

["]?, o;B (, , 1 ./ 2 )^Cr 2 -«- 2 "o||D„|| 2 2(B( ^ 3i . /4)) . 

By using the second Korn inequality, we get Caccioppoli's inequality for u, 
that is, for any A e R" and < p < r/2, we have 

f |Du| 2 dx<Cp- 2 (l + (diamQ) 2 ) f \u-A\ 2 dx. 

JB(x ,p) JB(x ,2p) 

Then, we get the estimate ( 12.231 from 1 12.251 . 



27 



Next, we prove the estimate ( 12.291 . We first establish a global version 
of 1161 Lemma 4.3]. For A > 0, we denote by L 2 ' A (U) the linear space of 
functions u 6 L 2 (U) such that 



IMta*(iQ = i SU P r A I |w| 2 dxl 

K=U JB(z,r)nU 



V0<J'<diam (J 



Lemma 6.1. Let f e L 2 ' A (B n Q)' ! , where A > and B = B(x ,R) with x e Q 
and < R < diam Q. Suppose u is a weak solution of 

[ Lu = f in BnQ, 
ii = on B n D, 
| t(w) = on B n N. 

If we assume H3, then, for < y < }>o ■= min(A + 4, n + 2/.( ), we «az?e 

r 2 "' f |Du| 2 <fa<c(f |Du| 2 rfx + ||/|| 2 ) (6.2) 

Js(i/)nQ Venn / 

uniformly for all x e \B fi Q and < r < R/2. Here, \B = B(x ,_R/2) and C is a 
constant depending only on n, k\, k 2 , j.io, A , A, y, and Q. We may take y = y in 
{6.2\ ify < n. Moreover, ify < n, then u e L 2 'i'(jB n Q) and 

W u Wl^ (ifinQ) ^ C(||u|| L 2 (Bnn) + ||Du|| L 2 (Bnn) + H/lb.A( B nQ)) ■ ( 6 - 3 ) 

Proof. Let x e B(x ,R/2) n <?Q and < r < r A (R/2), where r > is such 
that 

So) for some y e Q and 5 > 0. (6.4) 

xedn 

Denote 

N := Nn B(x,r), D := d(B(x,r) n Q) \ N 
and let z; be a unique weak solution of the problem 

Lv = f in B(x, r) n Q, s = on D, t(») = on N. 

By testing with v and using the Sobolev inequality, we get 

f |e(z;)| 2 dy<cf f \f\« dy) '(f |upiy) 

JB(i/)nu \Jfl(i,r)nQ / W8(.v/)nQ / 

where 1/p + l/q = 1 and q = 2n/(n + 2) if n > 3 and = 2/(a + 1) if n = 2, 
where a e [j.i , 1). We extend » to Q by setting v = in Q \ B(z, r). Note 
that v e W^Q)" and Dv = on Q \ B(x, r). By the Sobolev inequality and 
Holder's inequality, we then obtain 

lk(f)ll 2 2(Q) < Cr"l|/|| t 2 (Bnn) ||z7|| w i,2 (Q) . 

By the assumption i6A\ , we have 

c||D»|| L2(n) < ||£(w)|| L 2 (Q) , IMI w i, 2( n) < Q\Dv\\ L 2 (a) , 
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for some constants c, C > independent of v, x, and r. Therefore, we have 



f 

..'[!( 



\Dv\ 2 dy <Cr A+2a "-f" 2 



L 2 '' 1 (BnQ) ' 

S(i/)nQ 

Let w \— u — v. Then, w is a weak solution of Lw = in B(x, r) n O. For 
< p < r, we have 

( |Dn| 2 dy<2 f \Dv\ 2 dy + 2 \ \Dw\ 2 dy 

Js(x,p)nD JB(x,p)na JB(x,p)nn 



\ ' ' JB(i,r)nQ 



Thus, by (TTJ Lemma 2.1, p. 86], we obtain (62) for x e B(x ,R/2) n <?Q 
and < r < r A (R/2). The general case of the estimate \6.2\ is obtained by 
combining this case and the interior case, which is already covered by 1161 
Lemma 4.3]. The estimate J6.3> is an easy consequence of the estimate J6.2I 
and the Poincare's inequality; see [ll] Proposition 1.2, p. 68]. ■ 

Let u be a weak solution in Y 2 (Q n <3)" of <2.28t , where Q = Q_(X ,r) 
with X = (x , f ) e Q and < r < V^o A diam Q. Then, we have (see fl6l 
Lemma 4.2]) 



ess sup f \u t \ 2 (x,t)dx < Cr~ 6 J J |K| 2 dX, 

-(r/2) 2 <(<f JB(,xo,r/2)ra JjQnQ 
ess sup I |D x M| 2 (x,t)dx< Cr" 4 I I |w| 2 dX. 

-(r/2) 2 <f<t JB{x ,r/2)r,Q JjQnQ 



( -(r/2) 2 <t<( ( 

Also, we have (cf. (U Lemma 8.6]) 

I |u-A| 2 rfX< Cr 2 f |D x u| 2 dX; A := + udX. 

jQnQ jQr\Q J Qua 

By using Lemma[6j]and the above inequalities, we repeat the proof of tl6l 
Theorem 3.3] with obvious modifications to get 

[»W/%iana * A a r-^-t" +2 W 2 ||«||^. 2(QnQ) . (6.5) 

Finally, we show that the above estimate and the condition J2.26t implies 

|«| .i QnQ < ^ir-C +2 ' /2 ||«||^ 2(Qn Q). (6.6) 

For Y e j Q n Q and Zefl satisfying [Z - Y\& < r/4, we have 

| M(Y) |2 < 2 l-4 Wr 2, 1[M] 2 + 2 | M(Z) |2. 



By taking the average over Z and using {6.5) , we get 

| M(Y )p < [?-^A\ +p- 1 2 2 " +5 )r-"- 2 \\u\\% 2(QnQ) . 
By adopting a covering argument, we obtain \6.6l , and thus (Z29). 
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6.2 Proof of Lemma [5^91 

Recall that u is the weak solution of the problem \5.6i , where / is understood 
as an element of L 1 ((fl / b); V). We define 

f(x, t) := f(x, t) - n K [f(; t)](x), (x, t)eQ:=Qx (a, b). 

Note that / = J in L 1 ((fl, fc); V). Also, for all t e (a, fr), we have 

II^/(vOIIl»(Q) < Q\nnf(;t)\\ LHa) < C[l/(vf)lb(Q) < C|Q| 1/2 ||/||^ (Q) , 

where we used the fact that all norms in H are equivalent and that tcr is an 
orthogonal projection. Therefore, we have 

ll/II^CQ) < C||/||^ (Q) . (6.7) 
Let u be the weak solution in % {Of of the problem 

-ttf — Lu = f in Q x (a, b) 

u = on D x (a, V) 

t(u) = on N x (a, b) 

u{-, V) = on Q. 

Then it is easy to see that u(; t) e V for a.e. t e (a, b) and thus u is a weak 
solution of the problem i5.6i . Therefore, by the uniqueness, we conclude 
that u = u. In particular, u is a weak solution in ^(Q)" of —u t + Lu = f. 
Because of the inequality \6.7i , it is then enough to prove the second part 
of the lemma and we refer to |3] Section 3.2] for its proof. ■ 
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